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ABSTRACT 

Recent observations of galaxy clusters have shown that environmental effects apparently associated 
with the cluster begin to lower the star formation rates of galaxies at distances as great as three times 
the cluster virial radius. These observations may indicate preprocessing of cluster galaxies in groups 
or in the cluster core for galaxies on highly elliptical orbits, but may also imply that the environmental 
effects due to the cluster are directly affecting galaxies on their first infall. To explore these issues, we 
investigate different models of ram pressure stripping as it acts on satellite galaxies in clusters, and 
compare to observations of the radial star formation gradient in clusters. We calculate the location 
of the accretion shock around model clusters, and use this as the radius of onset of ram pressure 
stripping in the Galform semi-analytic model of galaxy formation. Comparison of the results of our 
model, and previously considered, simpler ram pressure models, with recent observations indicates 
that current data is unable to strongly discriminate between models of ram pressure stripping due to 
the complex interplay of preprocessing effects at work. However, future observations of a larger sample 
of clusters will likely be able to place stronger constraints on the process of ram pressure stripping 
and its role in shaping radial trends in and around clusters. 

Subject headings: cosmology: theory, galaxies: clusters: general, galaxies: evolution, galaxies: clusters: 
intergalactic medium 



1. INTRODUCTION 

It is widely known that the star formation rate (SFR) 
of galaxies depends on the density of their environment, 
and specifically that galaxies within galaxy clusters tend 
to be forming stars at a lower rate than comparable field 
galaxies. Since star formation is fueled by gas, this leads 
to the conclusion that galaxies in clusters tend to have 
less gas from which to form stars than their counterp arts 
in less dense environments ()Butcher fc Oemieilll978[ ). 

Several mechanisms have been proposed to explain the 
observed trend of star formation with environment. It 
is well known that spiral galaxies tend to be bluer, and 
therefore forming stars at a higher rate, than elliptical 
and SO galaxies, the so-called morphology-density rela- 
tion (pPressler 1980), and that spiral galaxies tend to be 
rarer in the centers of galaxy clust ers than early typ e 
galaxies (jWhitmore et al. 1993: Ba rkhouse et al.l 120091 ). 
It is clear, therefore, that morphological transformation, 
driven by interactions such as mergers or multiple weak 
gravitational e ncounters with oth er satellites (galaxy ha- 
rassment, e.g. iMoore et al.lll996D . likely plays a role in 
driving the observed quenching of star formation in clus- 
ters. However, this cannot completely explain the trend 
as it has been observed that even among galaxies with 
the same bulge to disk ratio, field galaxie s are forming 
stars at higher rates than cluster galaxies (jBalogh et al.l 
[1999). 

Gas may be removed from galaxies passing through the 
dense intergalactic medium of a cluster by means of ram 
pressure stripping (RPS). As originally envisaged, in this 
mechanism the interstellar medium (ISM) of a satellite 
galaxy is lost due to interaction with th e dense gas of 
the host halo through which it is moving ("Gunn fc Gott 
■1972; Abadi et al. 1999; Quilis et al. 2000; Kenne v et al. 
2009f ). Alternatively, a satellite moving through less 



dense material may lose just a portion of its diffuse atmo- 
sphere of hot halo ga s (colloquially known as starvation 
or str angulation, e.g. iBekki et al.l 1200^ iMcCarthv et al.l 
l2Q08f ). Finally, it is inevitable in a hierarchical Universe 
that many cluster galaxies were previously members of 
lower mass groups of galaxies, and environmental ef- 
fects such as strangulation in these lower density environ- 
ments may have begun the quenching of the star forma- 
tion of these galaxies before they we r e accreted onto the 
cluster (iZabludoff fc Mulchaevlll998l: iMcGee et al.l l2Q0l 
iKawata fc Mulchaevll2008l: lFuiitall200^ 

It has traditionally been assumed that these environ- 
mental effects only occurred once satellites had fallen 
through the virial radius of a host halo. However, re- 
cent observational studies of satellite galaxy SFR versus 
cluser-centric radius have found that the SFR remains 
depressed relative to the field out to 2 — 3 times the virial 
radiu s of clusters fBalogh et al.l 120001 : iVerdugo et al.l 
120081 iBraglia et al . 2009) . Three explanations have 
been suggested for the radial extent of environmental ef- 
fects: that group-scale effects in the locality of the clus- 
ter are having a large impac t on these cluster galaxies 
(IZabludoff fc Mulchaevl[l998l ). that many satellite galax- 
ies follow highly elliptical orbits which take them close 
to the cluster center, where they experience st rong RPS 
and then back out to large radius iBalogh et al 
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and that the radius at which cluster environmental ef- 
fects begins to become important is further out than the 
virial radius. 

In fact, it is already known that using the virial radius 
as the location at which environmental processes related 
to the cluster begin is only an approximation. In cur- 
rent structure formation scenarios, dark matter halos are 
formed through gravitational collapse when an overdense 
patch of dark matter stops expanding with the universe 
and collapses. The surrounding gas then falls into this 



potential well. If the scale of the halo is large enough so 
that the gravitational dynamical time of the halo is much 
less than the cooling time of the gas, then the gas will 
form an accretion shock, where the kinetic energy of the 
infalling gas is converted into thermal energy, heating the 
gas t o the virial temperature of the halo ( White fc Frenk 
I1991I) . These are the so-called "virial shocks," as they are 
expected to form near the virial radius of the halo. 

Accreti on shocks wer e predi cted in a cosmological sce- 
nario by iBertschinged (|1985| ). who found self-similar, 
spherically symmetric analytic solutions of collisional 
gas falling into a density perturbation in an Einstein- 
de Sitter universe. More recently, models of large- 
sca le, spherical accret i on sh ocks have been deve loped 
by iFurlanetto fc Loebl pOOl and|Bar^ii^ (|200l . and 
three-dimensional cosmological N-body and hydrody- 
namica l simulations have confir med the e xistence of such 
shocks (iKes het et al. 2003; Mol nar et aUWOO^ ). Further, 
iMolnar et al. (2009) has predicted that the extent of such 
shocks should be observable in the Sunyaev-Zel'dovich 
effect of clusters with next-generation radio telescopes 
such as ALMA. While virial shocks may be unstable, 
and so not survive, around low mass halos, they are ex- 
pected to be an inevitable result of s tructure formation 
for halos above a few t imes 10^^ Mq (|Birnboim fc Dekell 
[200aiKeres erall[2005l) . including groups and clusters of 
galaxies. When applied in the Galform semi-analytic 
struct ure formation code , the c riterion for shock stabil- 
ity of Birnboim fc DckeJ! ()2003[ ) gives the similar result 
that stable shock s form only in halos mo re massive than 
around IO^^Mq enson fc BowedlMoh . 

To determine the cosmological importance of accretion 
shocks, we apply this well-developed theory to galaxy for- 
mation, and in particular to the environmental effects on 
galaxies in clusters. Semi-analytic models of structure 
formation are a very powerful tool to investigate the im- 
pact of this additional physics on structure formation, as 
their relatively light computational requirements allow 
us to probe the efi^ects of accretion shocks on large sam- 
ples of galaxy clusters, therefore complementing the hy- 
drodynamical simulations mentioned above. Such semi- 
analytic models have proven to be very successful probes 
of galaxy formation; they have well re produced the lu- 
minosity function of galaxies locally ("C ole et all 120001 : 
iSomerville et al.ll2001t [B augh et al. 200$, and at high 
redshift fKauffman n et al.ll999h) . The addition of astro- 
physical effects such as supernova feedback, a photoioniz- 
ing background, and environmental effects have substan- 
tially im proved the fit at low luminosities ([Benson et al.l 
I2002allbl : ISomervilld [200^ . while the recent addition of 
feedback and heating effects to counter large-scale cool- 
ing flows, such as heat conduction and feedback from 
active-galactic nuclei, have produced faithful matches 
to the high-lumino s ity tail of the luminosity function 
(ICroton et al.l 120061: [Bower et al.l [20061: ISomerville et al.l 
l2008f ). The colors of satellite galaxies have also recently 
been ma tched by the add ition of a more detailed model 
of RPS (Font et al. '20081; see also iLanzoni et al.l [2005 
[Weinmann ct al. 2010). 

In this work, we consider the physical and observa- 
tional consequences of the onset of environmental ef- 
fects at the radius of the accretion shock rather than 
around the virial radius as has been previously assumed. 
To study this problem we apply RPS including accre- 



tion shocks to the clu sters in the Millennium Simulation 
([Springel et al.l[2005l) using the Galform semi-analytic 
model of structure formation. In section [2Tl we describe 
the Galform model that we use as the basis of our 
calculations, while in section [2.2! we briefly describe our 
calculation of the location of the accretion shock and 
its implementation in the semi-analytic model, giving a 
more detailed description of our calculation in the Ap- 
pendix. In section [3] we present the cluster galaxy prop- 
erties produced by our model, and compare them with 
other models of RPS. In section[4|we compare our model 
with recent observational results of the SFRs of cluster 
galaxies, and in section[S|we discuss the conclusions that 
we draw from these results for the observability of envi- 
ronmental effects in clusters. 

2. SIMULATIONS 
2.1. The GALFORM Semi-Analytic Model 

We use as the basis of our work the Galform semi- 
analytic model of galaxy formation as described by 
iBower et al.l ([20061 ). to which the reader is referred for 
a detailed de scription of t he m odel formu lation. As 
descr ibed in IBower et al.l ([20061 see also iHellv et all 
l2003f ). we apply this model to the detailed merger histo- 
ries of the entire vol ume of the Millennium Simulation 
JSpringel et"all I2005D . Our model uses the same cos- 
mological parameters as that simulation: fim = 0.25, 
Qb = 0.045, A = 0.75, and h = 0.73 at z = 0. 

As we are working with N-body merger trees, there 
are situations in w hich halos decrease in mass with time 
([Hellv et al.ll2003D . as a result of unbound particles incor- 
rectly being tagged as halo members for example. The 
original implementation of our semi-analytic model was 
not well equipped to deal with mass loss in halos, so 
these N-body merger trees were artiflcially forced to con- 
serve mass. However, in this work we utilize the merger 
tre es without req u iring mass conservation, as described 
by I Stringer et al.l ([20101 ) , as this is a fairer representa- 
tion of the true behavior of the N-body simulation and 
is important for this work in which we utilize halo mass 
growth rates to compute virial shock radii. We flnd that 
relaxing the requirement of mass conservation results in 
an increase in the mean stellar mass content of galaxies 
of around 0.3 dex. This chang e is larger than expected 
from the work of IHellv et alT ([20031 ) due to the sensi- 
tive nature of the AGN feedback included in our current 
mode l (but which was not present in that of IHellv et al.l 
[2OOI . 

We therefore found it necessary to adjust a single pa- 
rameter of our model relative to that of IBower et al.l 
(2006) to retain a good flt to the observed local bj- 
and K-band galaxy luminosity functions. Specifically 
we adjust the parameter acooi, described bv IBower et al.l 
([2006). which determines the halo mass scale above 
which AGN feedback becomes effective. We find that 
increasing acooi from 0.58 to 0.9 reduces the lower mass 
limit for AGN heating to become effective from a few 
times IO^^A/q in the Bower model to a few times IO^^Mq 
in this work, and therefore reduces the mean stellar mass 
of galaxies. This change tends to cancel out the change of 
stellar masses caused by not conserving mass in merger 
trees (as was done in IBower et al.|[2006D . and brings our 
model back into agreement with the local luminosity 



functions. 

The prescription for the treatment of reheated gas in 
sateUite ga l axies that we adopt is the same as that of 
I Font et aD (|2008f ). to which the reader is directed for 
complete description. In brief, gas in a satellite galaxy 
that is reheated by supernovae or AGN feedback is trans- 
ferred to the hot halo of the satellite, from which it is 
transferred to the host halo as the satellite is ram pres- 
sure stri pped. 

In the lBower et al.l (|2006( ) model, a galaxy that is iden- 
tified as part of a friends-of-friends group of a more mas- 
sive halo is considered its satellite. These tend to roughly 
correspond to satellites within the virial radius of the 
halo, although the prescription is by no means spheri- 
cally symmetric. The satellite has all of its hot gas in- 
stantaneously stripped away, leaving the cold gas in the 
ISM of the galaxy but removing the source of gas to be 
accreted onto the galaxy and form stars. Although this 
model has been quite successful in reproducing the lu- 
minosity function and star formation history of galaxies 
and their evolution, it fails to reproduce the colors of 
satellite galaxies, tending to predict them to be redder 
than is observed. 

The challenge to reproduc e the correc t color s of satel- 
lite galaxies was taken up by iFont et all (|2008D , who im- 
plemented a more nuanced approach to the RPS of satel- 
lite galaxi es bas e d on the hydrodynamic simulations of 
iMcCarth v et al.l ()2008[ ). Keeping the definition of satel- 
lite galaxies as those belonging to the friends-of-friends 
group of a more massive halo, in their model each satellite 
is assigned an orbit as s umin g the velocity distributions 
determined bv iBensonI ()2005[ ). and from this they calcu- 
late the maximum ram pressure exerted on the satellite 
by the host halo and galaxy, which occurs at the pericen- 
ter of its orbit. The radius at which this maximum ram 
pressure is equal to the gravitational restoring force per 
unit area of the satellite is termed the stripping radius, 
and all of the hot halo gas beyond this radius is stripped 
at the moment that the satellite crosses the virial radius 
of the host halo. This calculation was repeated (possibly 
resulting in more gas being removed from the satellite) 
every time its host halo doubled in mass since the previ- 
ous ram pressure calculation. This model provides a less 
extreme implementation of RPS, and manages to match 
the colors of satellite galaxies by allowing them to ac- 
crete hot gas and remain blue for a longer time after 
being accreted. 

However, the ram pressure model of lFont et all (pOOl 
continues to use the virial radius as the location around 
which a satellite begins to feel ram pressure from its host 
galaxy. In fact, ram pressure forces begin to be felt by 
the infalling galaxy when it passes through the accretion 
shock, at which the cluster gas temperature, density and 
pressure discontinuously increase. The accretion shock 
can be up to twice as far from the host galaxy as the 
virial radius, potentially significantly altering the effect 
of a massive host halo on nearb y galaxies. Additionally, 
the model of iFont et al.l (|2008| ) uses randomly assigned 
orbital parameters for satellites to compute their orbit 
and, therefore, the ram pressure force that they experi- 
ence. 

2.2. Implementation of Accretion Shocks 



In this work we calculate the location of the accretion 
shock of halos and use this as the radius at which the 
RPS of satellite galaxies begins, thereby more completely 
modeling the environmental effects of a host halo on its 
satellites. We calculate the radius of the accr etion shock 
with a model based on the calculations of iVoit et aTl 
(|2003l hereafter V03), with a few assumptions relaxed 
to obtain a more accurate accretion radius in a wider 
range of situations. We use the method of iFont et all 
(j2008| ) to implement the RPS of a galaxy once it comes 
within this radius of the cluster, but use the actual orbit 
of the satellite (taken directly from the N-body simula- 
tion) to compute the ram pressure force experienced at 
each timestep of our calculation. 

The method of accretion shock calculation of V03 is 
an approximate solution with many simplifying assump- 
tions. For example, it assumes smooth, spherical accre- 
tion, an assumption that is known to be quite incorrect 
in the context of hierarchical galaxy formation. For this 
reason, there is uncertainty as to the accuracy of the 
predictions of this model for the hydrostatic structure 
of cluster gas. For example, cosmological hydrodynamic 
simulations have shown that the hot gas halo of a clus- 
ter can extend well beyon d the location of the accretion 
shock ()Frenk et al.l 119991) . pointing to far more compli- 
cated physics than is included in the simple V03 model. 
This discrepancy is likely to be small, as the relative ve- 
locity between infalling satellite galaxies and accreting 
gas, which determines the ram pressure along with the 
density, will tend to be low outside the shock radius since 
they are both falling into the host halo and feel little 
pressure, while within the shock radius the accreted gas 
becomes nearly stationary creating a large relative ve- 
locity with the satellites. However, as it has been shown 
in numerical simulations of the hierarchical dark matter 
and gas evolution of galaxy clusters that the profiles of 
the outer regions of clusters tend to agree w ith smooth 
accretion models (|Borgani fc Kravtsovll2009[ ). we use this 
simplified model in this work. 

The V03 calculation of the accretion shock radius also 
assumes that the accretion shock is always perfectly 
strong, so that the Mach number of the shock approaches 
infinity. We relax the strong-shock assumption to allow 
the accretion shock to have any strength, which we find 
to be justified as even some massive clusters have accre- 
tion shocks with Mach numbers M ^ 10. We also correct 
for the neglect of an integration constant in the hydro- 
static profile derived for the clusters, which we found to 
be non-negligible. 

Briefly, the calculation of the accretion shock proceeds 
as follows: using the shock jump conditions and the as- 
sumption that the total accreted gas mass must be con- 
tained within the accretion shock, we obtain a simple 
hydrostatic model of the cluster gas. Using this model 
we derive an equation for the accretion radius in terms 
of the halo mass, mass accretion rate and halo concen- 
tration, as well as the Mach number of the accretion 
shock. By also assuming that the cluster gas has adiabat- 
ically contracted from a temperature of Tiqm = 3000K 
in the intergalactic medium, we simultaneously solve for 
the Mach number of the accretion shock and its position. 

The calculation on which our model is based is de- 
scribed in Appendix A of V03, and a more detailed de- 
scription of our generalized version of the calculation can 
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Fig. 1. — Stellar mass (M«, top), cold gas mass {Mgas, middle) 
and the gas fraction, Mgas /Ms, (bottom), averaged for satellite 
galaxies in radial bins, as a function of the cluster-centric radius 
divided by the virial radius of the cluster. The black dashed line 
shows the average position of the accretion shock in these clusters, 
while the blue, gre e n and red solid lin es show the results for the 
IBower et all J^TOl) . IFont et al.l 1 12003) . and Shocks models. The 
blue, green and red dotted lines show the Bower, Font and Shocks 
models without ram pressure stripping, respectively. The crosses at 
zero radius indicate the average properties of the central galaxies. 

be found in the appendix of this paper. 

3. CLUSTER GALAXY PROPERTIES WITH DIFFERENT 
MODELS OF RAM PRESSURE STRIPPING 

We first compare the properties of galaxies in our 
model, in which the stripping of gas begins at the ra- 
dius of the accre tion shock, with the similar model of 
I Font et all ()2008l ) in which these e ffects begins at the 
virial radius, and with the model of IBower et al.l ()2006f ) 
in which all of a satellite galaxy's hot gas, not just the 
maximum amount as determined by the parameters of 
its orbit within the host halo, is stripped away at the 
virial radius. We will call these models the Shocks, Font, 
and Bower models respectively. 

In Figure [T] we show the resulting stellar mass, cold gas 
mass and gas fraction profiles of cluster galaxies, exclud- 
ing the central galaxy, averaged within cluster-centric ra- 
dial bins. The values of the central galaxies are plotted, 
for comparison, as points at zero radius. Here we de- 
fine galaxy clusters as halos whose mass is greater than 
10^^h~^MQ, and we include all galaxies in each radial 
bin in the averages. It should be noted that the choice of 
which galaxies to include makes a large difference in the 
properties that are observed, as can be seen in the signif- 
icantly different trends obtained in figures [3] and |4] where 
we have selected a different set of galaxies to compare to 
observations. 

We compare our Shocks model (red solid lines), with 
the Font (green solid lines) and Bower (blue solid lines) 
models. It should be noted that there are other pa- 
rameter differences between these models: the Shocks 
model has had a parameter associated with AGN feed- 



back adjusted from its Bower model value to match the 
local galaxy b,j- and X-band luminosity functions, while 
the Font model has a different value of the metal yield 
to match the zero point colors of the red and blue se- 
quences. Thus, the difference between the Bower model 
and the Shocks and Font models is due both to their 
differing treatments of RPS, and to their different phys- 
ical parameters. To assess the relative contributions of 
each of these components, we also consider the Shocks, 
Font and Bower models without RPS, thus isolating the 
effects of the parameter changes (red, green and blue 
dotted lines) . We also show the position of the average 
accretion shock radius of these clusters as a black dashed 
vertical line, and the properties of the central galaxies as 
crosses. When these values are very different from those 
of the satellite galaxies, we leave them out of the plot as 
the satellite galaxy trends are more interesting for the 
purposes of this work. 

In the top panel of Figure [TJ we see that in all of the 
models there is a peak in galactic stellar mass in the cen- 
ters of clusters; as the central galaxy has been excluded, 
this shows an increase in stellar mass of the innermost 
satellites in all models. This is reasonable, as we ex- 
pect more massive satellite galaxies to sink deeper into 
the potential well of the cluster. Further, other than a 
slight rise in stellar mass towards the center of clusters in 
the Font and Shocks models, the average stellar mass of 
galaxies remains more or less constant with radius, with 
the Shocks model having the lowest average level of stel- 
lar mass, the Bower model slightly more, and the Font 
model having the most stellar mass of all. The flatness of 
these curves indicates that most of the stars of satellite 
galaxies were formed before they merged with the clus- 
ter. Note that the central galaxies in all of the models 
have more than an order of magnitude more stellar mass 
than the satellite galaxies, as would be expected since 
they are at the center of the potential well of the cluster 
and are not subject to RPS. 

The distribution of cold gas mass, in the middle panel 
of Figure [U shows a different trend. The Bower model, 
with its more extreme removal of gas in clusters, shows 
the least cold gas in satellite galaxies, with a very distinct 
drop in gas mass at around 1.5 times the virial radius. 
This drop in cold gas mass is quite steep, with its slope 
determined by the timescale of satellite orbits as com- 
pared to the timescale on which cold gas is made into 
stars in the satellite galaxies. Interestingly, the Shocks 
model shows a very similar trend with slightly more gas 
mass inside the radius of the accretion shock, as would 
be expected. That the Bower model also shows a drop 
at around this radius, which happens to be the aver- 
age accretion shock radius, is interesting, and points to 
the 'preprocessing' of satellite galaxies, in which satel- 
lite galaxies experience weaker environmental effects as 
members of smaller groups of galaxies before merging 
with the cluster. The Bower model is likely to exhibit 
stronger group effects than the Shocks and Font mod- 
els, since its RPS efficiency is always very high, while 
the RPS of the Shocks and Font models depend on the 
density of the halo intracluster medium (ICM). Thus, we 
see that these two very different models of environmental 
effects give a similar qualitative prediction for the radial 
dependence of satellite gas mass, indicating that if such 
a trend in gas mass is observed, we cannot distinguish 



between a model whose RPS starts at the accretion ra- 
dius and a much stronger RPS model which starts at the 
virial radius. 

The Font model predicts a higher level of gas mass at 
all radii, as is reasonable due to its less harsh RPS im- 
plementation, while as expected the models without ram 
pressure exhibit significantly higher gas mass at all radii 
than the other models, since it is RPS that is mainly 
responsible for the sharp drop in gas mass with decreas- 
ing cluster-centric radius. It is reasonable that there is a 
slight drop-off in gas mass towards the center of the clus- 
ter even in these models, since in these models the ac- 
cretion of new gas from the intergalactic medium (IGM) 
onto the hot gaseous halo of the satellites is suppressed 
though none of their halo gas is removed. We see that the 
central galaxies have higher gas mass than satellites, as 
they are not affected by RPS and also gain gas through 
mergers. The higher level of gas mass in central galax- 
ies in the models without RPS is likely due to merging 
satellites, which have not had gas removed by RPS and 
therefore give more gas mass to the central galaxy. 

Finally, in the bottom panel of Figure [U we can see 
the average ratio of cold gas to stellar mass in satellite 
galaxies. It is quite striking that, despite the differences 
in gas mass and stellar mass when averaged separately, 
the average of their ratio is nearly the same in all of 
these models, with the only significant difference occur- 
ring between the Bower model and the Font and Shocks 
models at less than half of the virial radius. The general 
behavior of this quantity, in all of the models, is a grad- 
ual increase in gas fraction moving into the cluster until 
the virial radius, at which the gas mass declines precip- 
itously. This peak is likely due to the biased sample of 
galaxies located close to massive clusters. As they are in 
general more massive, they may well still be cooling gas 
and thus forming stars at rates higher than further away 
from the cluster. 

Further inwards, the Shocks and Font models reach a 
higher, roughly fiat, central gas fraction than the Bower 
model, as expected due to their less extreme removal of 
gas in cluster members. We also see that the gas fraction 
reaches only a slightly larger value in the runs without 
any RPS, indicating that at least the subtle ram pressure 
of the Shocks and Font models makes little difference in 
the gas fraction of cluster galaxies, while the implemen- 
tation of a Bower-type complete removal of gas makes a 
far larger difference. Note, finally, that the spike in the 
Bower model gas fraction very close to the center is due 
to a single galaxy in our sample, which formed in the 
final timestep and so exhibits an unusually low stellar 
mass and therefore a very high gas fraction. Its loca- 
tion deep within the cluster is surprising given its very 
recent formation and it may represent a flaw in the halo 
detection or tree building algorithms. We leave it in the 
sample in any case for completeness. Finally, the gas 
fractions of central galaxies in these models are all well 
below those of satellites, which is reasonable as cluster 
central galaxies tend to be ellipticals which have little 
cold gas although they may have large hot gas halos. 

In Figure [21 we show observationally measurable quan- 
tities related to the implementation of RPS: the B— V 
color, SFR, and specific SFR of satellite galaxies in clus- 
ters. Each of these quantities, as in the previous figure, 
is calculated for each galaxy and then averaged within 
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Fig. 2. — B-V colors (top), star formation rate (middle), and spe- 
cific star formation rate (bottom)), averaged for satellite galaxies 
in radial bins, as a function of the cluster-centric radius divided by 
the virial radius of the cluster. The lines and models are the same 
as in figure [T] 

cluster-centric radial bins. 

In the top panel of Figure [2] we can see that in all 
of the models the satellite galaxies become increasingly 
redder moving towards the center of the cluster as ex- 
pected, with this reddening beginning around the virial 
radius and having the same general shape in all mod- 
els. We see that the Shocks model has bluer colors than 
the other models, and that the Font model has slightly 
redder colors than the Bower model. The Font model, 
with the most gas, has the most dust extinction, and it is 
this effect which makes it redder than the Bower model. 
As the separation between the Shocks and Font models 
is preserved when RPS is turned off, we see that these 
average color differences result from the star formation 
and AGN feedback parameter adjustments which were 
required to bring the models into agreement with the lo- 
cal luminosity function. As expected, without RPS we 
see much bluer colors in all of the models. The central 
galaxies are all much redder than their satellites, which 
follows from our knowledge that these galaxies tend to 
be large elliptical galaxies. 

In the middle panel of Figure [U we see the SFR varia- 
tion with cluster-centric radius. The SFR versus radius 
for the Shocks and Font models decreases moderately in- 
side the virial radius, in a similar manner to the gas mass, 
though not quite as steep. We see a distinct decrease in 
star formation in the Bower model towards the center 
of the cluster, due to its harsher implementation of RPS 
and thus lower gas mass. The Font model has more star 
formation at all radii with respect to the other models, 
while the Bower model always has the lowest SFR. This 
is reasonable as it mirrors the order of strength of RPS 
in these models. Also as expected, the models without 
RPS show higher rates of star formation than the other 
models. The central galaxies in these models have much 
larger SFRs than their satellite galaxies, which is reason- 



able given their larger gas mass. 

In the lower panel of Figure [2l we see the average spe- 
cific SFR, that is the average of the ratio of the SFR and 
the total stellar mass of satellite galaxies. The trends in 
this plot are very similar to those of the gas fraction plot- 
ted in the lower panel of figure [TJ This is no coincidence; 
in fact the SFR in Galform is determined such that 
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where e is the star formation efficiency, Tdisk is the disk 
timescale, and R is the fraction of mass going into stars 
that is recycled back into the interstellar medium. As 
with the gas fraction, a single outlier galaxy in the Bower 
model is causing the spike at low radius, and can be ig- 
nored. We see that the central galaxies in all models have 
significantly higher specific SFR, as expected. Also, the 
specific SFR is directly proportional to the gas fraction, 
and as e and R are constants and Tdisk does not depend 
strongly on environment, we see similar behavior in the 
specific SFR as a function of cluster-centric radius as we 
saw in the gas fraction. As with the gas fraction, it is 
very similar in all models and shows a very slightly rising 
level moving in towards the virial radius, at which there 
is a peak and then a sharp decline towards the center of 
the cluster. As with the gas fraction, the models with- 
out RPS reach a higher central value of the specific SFR, 
the Font and Shocks models are very similar, and the 
Bower model shows a much lower central specific SFR. 
The central galaxies in these clusters all show a much 
lower specific SFR than the satellites; this is directly re- 
lated to their low gas fraction, which is understood as 
these galaxies tend to be large ellipticals. 

4. COMPARISON WITH OBSERVATIONS 

We compare the results of our Galform model in- 
cluding accretion shocks an d the more extreme RPS 
model of iBower et al.l ()2006l) to recent observations and 
simulations. In section | 4 .1 j we compare with the 
simulations of 'Balogh et al.l ( 200C1[) and their observa- 
tional dataset selected from the CNOCl cluster red- 
shift survey (|Yee et all Il996( ) . In section 112] we com- 
pare the Shocks mo del to the cluster spectroscopy of 
iVerdugo et al.l (|2008l ). In both papers, cosmological pa- 
rameters A = 0.7, J7o = 0.3, and h = 0.7 are used. 
Despite the fact that Galform uses a different set of 
cosmological parameters as listed in section 12.11 we con- 
duct our comparison analyses using the same parameters 
as the observations so as to better reproduce their anal- 
ysis. In both comparisons, we plot the properties only of 
satellites in halos with masses greater than lO^'^ h~^MQ 
to more accurately mimic the selection of massive clus- 
ters in these two samples. 



4.1. Comparison to CNO Cl Cluster R e dshift Survey as 
Selected bv \Baloah etall \200(\ ) 

The cluster redshift survey CNOCl (lYee et al.l[T99fi[) 
provides spectra for member galaxies of 15 X-ray lumi- 
nous clusters. We compare our simulation s to the sam- 
le of twelve of these clusters presented in iBalogh et al.l 
20001) ■ which were selected to lie within a redshift range 
0.19 < z < 0.45 and have well-defined cluster cen- 
ters. Cluster members were selected based on velocity 
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Fig. 3. — Co mparison of the cluste r redshift survey CNOCl data 
as obtained by IBalogh et al] I I2000I ). with the Millennium results 
using three different implementations of RPS. The CNOCl data 
is shown as black squares a nd solid lines , while the extreme ram 
pressure stripping model of IBower et al.l lj2006) is shown in blue 
dotted lines, the more nuanced accretion shocks model of this paper 
is shown in red dashed lines, and this same Shocks model without 
ram pressure stripping is shown in purple dot-dashed lines. The 
vertical lines show the average virial radius (cyan long-dashed) and 
accretion shock radius (green short dot-dashed) for the clusters 
considered in our semi-analytic model. Note that the x-axis of 
this plot is the projected cluster-centric radius divided by r2oo as 
defined in the text. 

and magnitude cuts. The cluster-centric radii of galax- 
ies were measured relative to r200) the radius at which 
the interior density is 200 times the critical density. See 
IBalogh et al.l (|2000[ ) for a complete description of the se- 
lection criteria. 

To ensure accurate comparison to the data, we mimic 
th e observationa l tech niques as described in section 2 
of IBalogh et al.l (|2000f ) in our analysis of the Millen- 
nium/GALFORM models. We analyze clusters at red- 
shift z = 0.3, and cluster members are selected based 
on projected radial position, magnitude, and velocity. 
Specifically, we select galaxies whose projected cluster- 
centric radius is less than 2 r2oo , where r2oo is calculated 
assuming an NFW-type density profile, with a Gunn r- 
band magnitude greater than —18.8 -I- 51og/i at z= 0.3, 
and with velocity within 3 u of the average velocity of 
all previously selected galaxies, including the effect of 
Hubble expansion. The average virial radius (as defined 
by Galform using the overdensity of a spherical top- 
hat collapse model for this cosmology) of t he clusters we 
"obser ve" in this manner is 1.37 r2oo- As in IBalogh et al.l 
(l2000l) . we remove the central cluster galaxies from the 
sample. 

We calculate the SFR as was done with the observa- 
tions, by using the equivale nt width of the [O II] line 
as a direct indicator of SFR. IBalogh et al.l (|2000l ) find a 
prescription to determine SFR from observed line widths 
and luminosities that matches the relationship between 
these quantities that they see in their semi-analytic sim- 
ulations. In their prescription, the SFR is the product 



of this equivalent width with the rest-frame B-band hi- 
minosity, an extinction factor, and a normahzation con- 
stant chosen by comparison with their simulations. In 
Figure [3] we show t he comparison of the CNOCl data of 
iBalogh et al.l ([2000) with the results of taking our semi- 
analytically determined cluster galaxy luminosities and 
line widths and plugging them into this prescription. We 
note that this prescription tends to under-predict the 
SFRs of our galaxies as compared to the rates directly 
obtained from Galform, and that this effect is strongest 
for models with less RPS, more galactic gas and dust 
and thus more extinction. Thus, the SFRs inferred from 
this method agree fairly well with the directly computed 
rates for the Bower model, but are much lower than the 
directly computed rates of the Shocks model. This shows 
that prescriptions to determine SFR from observed quan- 
tities are quite model dependent, and in this case disagree 
with the relationship we see in Galform. 

The resulting curves are compared with the CNOCl 
data in Figure [S] Her e we plot three mo dels of ram- 
pressure stripping: the iBower et al.l ()2006f ) model, our 
new Shocks model, and this same Shocks model with 
the RPS turned off (blue dotted, red dashed and purple 
long-dot-dashed lines, respectively). We also show, with 
cyan long-dashed and green short-dot-dashed lines, the 
average virial and accretion shock radii of the clusters 
considered. The radius plotted on the x-axis is the pro- 
jected cluster-centric radius. Note that the SFRs shown 
in figure [3] are different from those in figure M due to 
their differing definitions of radius, different selection of 
galaxies, the different units and the linear scale in figure 
[3] as compared to the logarithmic scale in figure [H When 
taking these into account, the SFRs in the two plots are 
similar. 

We can see that all three models predict similar radial 
trends in SFR. The Bower model, with more extreme 
RPS, is quite close to the data in the innermost regions, 
but has a much steeper rise to larger radii than is ob- 
served, thus overproducing the SFR in the outer regions. 
The Shocks model rises more gradually in the outer re- 
gions of the cluster, bringing both of these models into 
better agreement with the data than the Bower model 
at large radii. However, we can see that, given the large 
error bars on the data, neither the Shocks model (with 
RPS) or the Bower model is a significantly better fit to 
the observations. This interesting result implies that, 
even if only RPS is acting to cause these gradients in 
SFR, from data such as these alone we cannot distinguish 
the details of the active RPS mechanism. Nevertheless, 
our Shocks model is marginally the best match to the 
data. 

4.2. Comparison to I Verduao et ai\ i2003i ) Cluster 
Spectroscopy 

We also compare with the sam ple of six galaxy clus - 
ters spectroscopically observed bv lVerdugo et al.l (|2008[ ). 
The star formation activity of cluster galaxies as a func- 
tion of cluster-centric radius was investigated using the 
equivalent widths of the Ha and [O II] emission lines, 
and a significant drop in star formation in cluster galax- 
ies relative to the field was observed. Fo r the details of 
the ob servations, the reader is referred to lVerdugo et al.l 
(|2008h . 

To mimic the observational selection and analysis, we 




Fig. 4. — Comparison of the I Verdugo et aLl I I2008I) cluster spec- 
troscopy with our semi-analytic results using three different imple- 
mentations of ram pressure stripping. The x-axis is the projected 
cluster-centric radius divided by r200, which is defined differently 
than in the previous section, as explained in the text. The solid 
black lines with dark grey shading show the observed cluster re- 
sults and their errors, while the solid grey lines surrounded by light 
gray show the observed results and errors for field galaxies. The 
models shown here are the same as in Figure [S] 

determine the mean cluster redshift and the cluster ve- 
locity disper sion using the bi-w eight location and scale 
estimators of' Beers et al .' (1990), and select cluster mem- 
bers as galaxies whose redshift is within 3cr of the mean 
cluster redshift. We use the projected cluster radius nor- 
malized to r2oo, where r2oo is here defined to be 



We further remove from our cluster galaxy sample 
those galaxies whose apparent I-band magnitude is above 
the survey spectroscopic limit / = 19.5, those whose line 
of sight velocity relative to the cluster center places them 
clearly outside of the cluster, and those whose absolute l- 
band magnitude is greater than the limit of Mj = —21.4, 
which was imposed on the observed data to treat clusters 
at different redshifts with the same luminosity restric- 
tions. 

The com parison of the observational results of 
IVerdugo et al. (2008) with our ram pressure implementa- 
tions is shown in Figured The solid black lines with dark 
grey shading show the observed cluster results and their 
errors, while the solid grey lines surrounded by light grey 
show the observed results and errors for field galaxies. To 
this we compare our models; as in Figure [3l we show the 
Bower model (blue dotted), and the Shocks model with 
(red dashed) and without (purple dot-dashed) RPS. We 
also show the average virial radius (cyan long-dashed) 
and accretion radius (green short dot-dashed) for the 
Millennium clusters considered. 

In Figure |4] we can see that, as in Figure |3l the average 
virial radius of these clusters is around 1.4 r2oo, while 



the average accretion shock radius is at 2.5 r2oo- This is 
just at the edge of the radial extent of the observations, 
and indicates that if our model of the accretion shock 
is correct, then the current observations are not probing 
the cluster beyond the radius of environmental effects. 

The top panel of Figure |4] shows the radial dependence 
of the [O II] equivalent width. We see that the three 
models predict very similar, nearly straight-line radial 
profiles, all of which are plausible given the size of the 
observational error bars. Therefore, these cluster galaxy 
[O II] data cannot reliably determine whether any form 
of RPS is causing the observed decline in SFR. However, 
all of the models tend to predict a shallower slope than 
is observed. 

In the lower panel of Figured we see the radial profile 
of the Ha line equivalent width. In this case, both the 
Bower and Shocks models are plausible with the given 
error bars, while the Bower model predicts a lower overall 
equivalent width and fits the data slightl y better. 

It c an be seen that the data from both iVerdugo et all 
(|2008D and lBalogh et al.l (|2000D have a local star forma- 
tion peak at around 0.5 r2oo- The statistical significance 
of these peaks, due to the size of the errors, is ques- 
tionable, and their correlation can not be immediately 
understood to be due to the same phenomenon, as r2oo 
was defined differently in these two papers. However, if 
these peaks do in fact represent a general feature of clus- 
ter members' star formation rates, then this effect is the 
product of physics which has not been included in the 
current model. 

5. DISCUSSION AND CONCLUSIONS: IMPLICATIONS FOR 
OBSERVATIONS OF RAM PRESSURE STRIPPING 

We have addressed the question of the nature of the 
environmental effects felt by satellite galaxies in galaxy 
clusters by implementing three different prescriptions of 
RPS onto the Millennium Simulation, using the semi- 
analytic galaxy formation model of Galform. The 
three model s imp lemented are the complete RPS of 
I Bower et al.l ()2006f ). in which all of the hot gaseous halo 
of a galaxy is removed when it is first identified as a mem- 
ber of the friends-of-friends group of a ni ore massive halo, 
the more nuanced model of iFont et al.l ((2008), in which 
hot halo gas is only stripped up to a maximum stripping 
radius determined by the orbit of the satellite, and our 
new model which incorporates the nuanced Font RPS 
model but differs from their model by using the radius 
of the cluster accretion shock as the location at which 
these effects begin and by utilizing the actual satellite 
orbit measured directly from the N-body simulation. 

Considering the results of applying these models to 
the detailed merger histories provided by the Millennium 
Simulation, we see that, in general, all of the models show 
similar trends in stellar and gas mass with cluster-centric 
radius, with differences in normalization due to the dif- 
ferent ram pressure models and due to the different pa- 
rameters adopted in the models to obtain a close fit to the 
properties of the local galaxy population. Further, the 
gas fraction and therefore the specific SFR is very simi- 
lar in all of the models considered, with differences only 
apparent in the regions less than half of a virial radius 
from the center of the cluster. We see that, in particular, 
the Shocks and Bower models predict very similar radial 
average gas mass profiles, with a sharp down-turn at 1.5 



virial radii, the average radius of the accretion shock in 
these clusters. This indicates that if we were to observe 
the gas mass trend in cluster galaxies, we would not be 
able to distinguish these two very different implementa- 
tions of RPS. 

As expected, we find that the B-V colors of satellite 
galaxies in all of the models become redder towards the 
center of the cluster, and that in all of the models we con- 
sidered these effects begin just outside the virial radius 
of the cluster. The onset of this reddening well outside 
of the virial radius, as well as the decrease in gas mass 
beyond than the virial radius, indicate either the prepro- 
cessing of satellite galaxies before their accretion onto 
clusters or the presence of satellites on highly elliptical 
orbits which have already passed through the central re- 
gions of the cluster. 

We have compared the result s of our accretion shock 
model of RPS and that of .Bower et al.l (12006.) with the 
ob servationa l resul ts of IVerdugo et al.l (I2008D and thos e 
of lYee et al.l (|1996D as selected bv iBalogh et ail (|2000D . 
We see in these comparisons that both of these mod- 
els are consistent with the observations, given their un- 
certainties. Thus, we conclude that current observa- 
tional data on the radial SFR gradient in clusters do not 
strongly discriminate between different models of cluster 
environmental effects. However, our models clearly in- 
dicate that the presence of RPS has a strong effect on 
radial trends within clusters and we expect that it will 
be possible to constrain the details of these models given 
future observations with larger samples of cluster galax- 
ies. 

Recent work on the semi-analytical modeling of satel- 
lite galaxy stripp i ng eff ects in clusters was carried out by 
iWeinmann et al.l (|2010l ) . As in this work, they have mod- 
ified previous models, in which all halo gas was stripped 
as soon as a galaxy becam e a satellite in a larger halo 
(|De Lucia fc Blaiz"otll2007D . to include a more nuanced 
form of gas stripping. They model cluster-environment 
effects entirely with tidal stripping, and presume that 
gas is lost from satellite galaxies in proportion to the 
loss of dark matter. The baryonic physics of RP S is ne- 
glected in the model of IWeinmann et aD ()201Ql) . while 
in the model of this work we neglect tidal stripping. Of 
course, in cosmological structure formation both of these 
mechanisms have an effect on the stripping of gas from 
satellite galaxies. However, using the virial scaling rela- 
tions, we find that the ratio of the forces of ram pressure 
and tidal stripping is given by 
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where A is a constant of order unity, so the force of RPS 
larger than that of tidal stripping, for a satellite at the 
virial radius, by a factor proportional to the ratio of the 
mass of the host halo to the mass of the satellite to the 
two-thirds power. By this argument then, at least at the 
virial radius RPS exerts a larger force on gas in a satellite 
halo, and so we can be confident that we have included 
the domi nant physical effect in th e present work. The 
results of IWeinmann et al.l ()2010l ) are compatible with 
our results, and they also find that several of their models 
are plausible given the error bars of the data with which 



they compare. 

In conclusion radial trends in galaxy properties around 
clusters can now be accurately predicted by the Gal- 
FORM model of galaxy formation. Measurements of these 
trends therefore have the potential to place strong con- 
straints on the processes of mass accretion and star for- 
mation, both of which are key components of our pic- 
ture of galaxy formation. We have not discussed the 
distribution of cluster gas in detail in this work, but it is 
clearly a key ingredient in any model invokin g RPS as a 
driver of cluster galaxy evolution. Recently, Bower et al.l 
dSlOS) described a more advanced calculation of cluster 
gas physics within Galform which aimed to match the 
X-ray properties of clusters. Future work in this sub- 
ject should clearly explore both cluster galaxy and X-ray 
properties in tandem to ensure realistic modeling of the 



cluster physics. This, coupled with larger samples of clus- 
ter galaxies would greatly improve the statistical power 
of this method as an important constraint on galaxy for- 
mation physics. 
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APPENDIX 

APPENDIX A: ACCRETION SHOCK CALCULATION 

The accretion shock model of lVoit et al.l ()2003[ ) upon which we base our implementation makes several assumptions. 
The cluster is assumed to be spherically symmetric and in hydrodynamic equilibrium, with an effective equation of 
state P{r) cx [p{r)]'''''' , where 7eff = 1.2, and with the actual equation of state of a free monatomi c gas, P = K(s)p^^^ , 
where K is a function of the specific entropy. The cluster potential is assumed to be of the NFW (jNavarro 
type, and the post-shock velocity of the gas is assumed to be negligible. 



With these assumptions, they obtain the following hydrostatic model of the gas within the accretion shock: 



T{x) = TAg{x) 
p(x)=p,[5(x)]i/(7c«-i) 
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Here, x = r jrmr is the ratio of the cluster-centric radius to the virial radius, pg and g\ are constants of integration, 
which can be specified by the constraint that the baryonic mass inside the shock radius is equal to the total baryonic 
mass of the halo /b M and the shock jump condition relating post-shock temperature to the incoming velocity, 
respectively. 

In our implementation of this model we relax two particularly unjustified assumptions, that the accretion shock is 
always strong, and that the constant g\ in equation IA4l is negligible. To relax these assumptions we apply the further 
(better founded) assumption that the gas accreting onto the halo is compressed and heated adiabatically as it makes its 
way to the shock radius from the mean density IGM. Using this assumption, and also assuming that the temperature 
of the mean IGM is 3000 K, we can calculate the temperature just outside the shock, from which we can obtain the 
Mach number of the shock M. . 

With this prescription for finding the strength of the sh ock M., where — ) ■ oo is the strong-shock limit initially 
imposed bv lVoit et al.l ()2003[ ). we use the shock relations of iLandau fc Lifshitd ()1959() to obtain an equation for gy. 
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Here, the ratio of the accretion shock radi us to the v irial radius is shown as Xac- We finally find that the equation 
for the accretion shock, equation A12 in Voi t et al.l (|2003i) . is modified to 
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where C is the shock compression factor, the ratio of the densities internal and external to the accretion shock. 
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and we use the variable p which is defined in iVoit et aD (|2003f ) to be 

. 4 2 

There is a further subtlety in the calculation of the mass accretion rate onto halos, which has a strong effect on the 
final accretion shock radius through its appearance in equation lAlGI We have explored several possibilities, but in 
this work we adopt a spherically symmetrized model averaged over the dynamical time of the halo. More specifically, 
we calculate the average mass accretion rate over the past one dynamical time of the halo, including all of the mass 
added to the most massive progenitor of the current halo through mergers or accretion. Since we use this mass as 
an estimate of the spherically accreted mass of the halo, we are making the assumption that all of this added mass, 
even that which was added through mergers, will distribute itself and cause accretion shocks much as the same mass 
of spherically symmetric accretion would. 



